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ABSTRACT 

Context. Charge transfer (or exchange) reactions between hydrogen atoms and protons in collisionless shocks of supernova remnants 
(SNRs) are a natural way of producing broad Balmer, Lyman and other lines of hydrogen. 

Aims. We wish to quantify the importance of shock-induced, non-thermal hydrogen emission from SNRs in young galaxies. 
Methods. We present a method to estimate the luminosity of broad (~ 1000 km s ) Lya, Ly8, Lyy, H/? and Pa lines, as well as the 
broad and narrow luminosities of the two-photon (2y) continuum, from existing measurements of the Ha flux. We consider cases 
of p = 0. 1 and 1, where p = T e JT p is the ratio of electron to proton temperatures. We examine a modest sample of 8 proximate, 
Balmer-dominated SNRs from our Galaxy and the Large Magellanic Cloud. The expected broad Lya luminosity per object is at most 
~ 10 36 erg s -1 . The 2y continuum luminosities are comparable to the broad Ho- and Lya ones. We restrict our analysis to homogenous 
and static media. 

Results. Differences in the Lya/Ha and LyS/Ha luminosity ratios between the /3 = 0.1 and 1 cases are factors ~ 2 for shock velocities 
1000 < v s < 4000 km s" 1 , thereby providing a direct and unique way to measure /?. In principle, broad, "non-radiative" Lya from 
SNRs in young galaxies can be directly observed in the optical range of wavelengths. However, by taking into consideration the 
different rates between core collapse and thermonuclear supernovae, as well as the duration we expect to observe such Lya emission 
from SNRs, we expect their contribution to the total Lya luminosity from z ~ 3 to 5 galaxies to be negligibly small (~ 0.001%), 
compared to the radiative shock mechanism described by Shull & Silk (1979). Though broad, non-thermal Lya emission has never 
been observed, these photons are produced in SNRs and hence the non-radiative Lyo- luminosity is a part of the intrinsic Lya spectrum 
of young galaxies. 
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1. Introduction 

Observations of galaxies at high redshifts have revealed a broad class of Lya-emitting galaxies at z ~ 3 to 5 (e.g., Tapken et al. 
2007). The Lya emission from these objects is reaching us as light in the visible spectral band, enabling their study using large, 
ground-based optical telescopes, which in turn permits detailed spectroscopic studies of these galaxies. Observations of quasars 
at z ~ 6 (e.g., Fan et al. 2006) have revealed heavy elemental abundances exceeding solar values. We know that at least some of 
the galaxies at z ~ 3 to 5 have high abundances of heavy elements, facilitating the formation of dust. In homogenous and static 
media, the dust particles impede the escape of Lya emission from gas-rich galaxies, due to the small mean free paths of the photons, 
low temperatures of the gas and ultimately high probabilities of absorption. In clumpy media, dust can enhance the escape of Lya 
photons relative to the continuum (Neufeld 1991; Hansen & Oh 2006). Broadening of Lya lines due to multiple scatterings is a slow 
process requiring a long diffusion time (though velocity fields in the interstellar medium may broaden the Lya lines and reduce the 
diffusion time). Hence, there is special interest in the physical processes that are able to naturally produce extremely broad wings in 
Lya lines, which may permit the photons to leave the host galaxy without requiring many scatterings (but see ^5). 

Among obvious mechanisms is the one at work in the unique massive binary SS433 (for a recent review, see Fabrika [2004]), 
with strongly blue- and redshifted Ha and H/3 lines, due to cooling and recombination of hydrogen in the baryon-dominated, 
precessing jet moving with velocity ~ 0.26c. Such objects are very rare — SS433 is the only such example in our Galaxy. More 
well-known Galactic sources of Ha emission with broad line wings are the supernova remnants (SNRs) of Type la, emitting due to 
charge transfer (or "charge exchange") reactions between hydrogen atoms and protons in the blast wave penetrating the low-density 
(~ 1 cm 3 ), ambient gas. The widths of the Ha lines correspond to Doppler broadening with velocities up to ~ 5000 km s _1 . The 
same process should produce not only Ha emission, but photons in the Lyman series of hydrogen as well. Recently, some of these 
SNRs were observed in Ly/3 using the FUSE spacecraft (Korreck et al. 2004; Ghavamian et al. 2007, hereafter G07). 

Knowledge of the cross sections of charge transfers to excited levels and excitation of the fast-moving hydrogen atoms permit us 
to find simple formulae relating the luminosities of SNRs in the broad Ha and Lya lines. The Lya line should have a similar spectral 
distribution to the observed Ha one in the broad wings, because the optical depth of the SNR for broad photons is negligibly small 
and the optical depth for coherent scattering (in the distant Lorentzian wings) in interstellar gas is low. 
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We compile the existing data for core collapse and thermonuclear SNRs, including SNR 1987 A (where the reverse shock is 
bright in the broad Ha line), and present their theoretically expected, broad Lya and Lyp 1 luminosities. For two objects, we present 
their expected broad Lyy, Hp 1 and Pa luminosities. Taking into account the supernova (SN) rates, the luminosities of the SNRs in 
Ha and the duration of their active phase (for the charge transfer mechanism described), we find that — even without discussing the 
cosmological evolution of the SN rates — the expected broad Lya is several orders of magnitude lower than the estimate of Shull & 
Silk (1979), who treated fully radiative SNRs with low metallicities and velocities (20 to 120 km s ). We come to the conclusion 
that the contribution of both core collapse and thermonuclear SNRs to the Lya luminosity of young galaxies is negligibly small. 
In $2l we gather a modest sample of 8 Galactic and Large Magellanic Cloud (LMC) remnants, and use them as a template for 
estimating the expected Lya, Ly/3, Lyy, Hfi and Pa production. In §0 we compute the Lya/Ha, Lya/Lyp 1 , LypTHa, Lyy/Ha, H/3/Ha 
and Pa/Ha luminosity ratios. We present our results in $4]and discuss their implications in $5] 



2. Galactic & LMC Remnants 

SNRs are the result of the interaction of SN ejecta with ambient matter. The nature of the interaction can be approximately cat- 
egorized into several stages (Truelove & McKee 1999, hereafter TM99; and references therein): the ejecta-dominated (ED) or 
freely-streaming stage; the Sedov-Taylor (ST) or self-similar stage; the pressure-driven snowplow (PDS) stage; and a possible, 
momentum-conserving snowplow stage (Cioffi, McKee & Bertschinger 1988). Many of the well-studied, young SNRs like Kepler, 
Tycho and SN 1006 are intermediate between the ED and ST stages; this has been corroborated by the numerical studies of TM99, 
who showed that there is no sharp transition between the two stages. The transition from the ED to ST stage occurs on a timescale 
tso ~ t C h', the characteristic timescale is 

5/677-1/2 -1/3 ,„ 
fch = fch.o mj. E 51 ' n , (1) 

where M e j = yii^Mq is the mass of the ejecta, E = £51 10 51 erg is the energy of the supernova explosion, and no is the density of 
the ambient medium (in cm -3 ). The coefficient in equation (Q~|l is f^o = 423 yrs (TM99). If one makes the argument that M e j (with 
density p = 1.4#ih«o) of mass is swept up in a time f c h, one instead gets f^o =186 yrs. The PDS stage occurs at 

fpos ~ 30f ch m e . E 5 \ n Q Q (2) 

after the explosion (Cioffi, McKee & Bertschinger 1988; TM99), where ( m is a dimensionless metallicity correction factor. More 
precise estimates for ?sd and ?pds are dependent upon the spatial density distributions of both the ejecta and the ambient matter. 

In the ED and SD stages, the emission from some SNRs is "non-radiative", meaning the timescale for thermal, radiative losses 
from the interacting gases is much longer than t c h. When the blast wave of the SNR slams into ambient gas consisting predominantly 
of hydrogen atoms, it emits in Balmer and Lyman lines consisting of abroad (~ 1000 km s -1 ) and a narrow (~ 10 km s _1 ) component 
(Chevalier & Raymond 1978; Bychkov & Lebedev 1979; Chevalier, Kirshner & Raymond 1980; Heng & McCray 2007, hereafter 
HM07; Heng et al. 2007, hereafter H07; G07; and references therein). These objects are known as "Balmer-dominated" SNRs. 
Positive detections of the line components are so far only from Galactic and LMC SNRs. Even though narrow Lya emission is 
produced, it is not seen due to interstellar absorption; broad Lya should be observecfl Non-thermal Ha and Lya emission has not 
been observed in studies of local starburst galaxies (e.g., Kunth et al. 2003). 

The narrow Balmer and Lyman lines are produced when the fast-moving ejecta directly excite stationary hydrogen atoms in the 
surrounding material. The broad lines are produced when the post-shock protons and atoms engage in charge transfer reactions, 
creating a population of post-shock atoms in broad velocity distributions known as "broad neutrals" (HM07; H07). In the frame of 
the observer, these broad neutrals move at a velocity vb S 3v,/4, where v s is the shock velocity (of the blast wave). For v s > 500 km 
s _1 , the broad neutrals can produce Lya that is blue- or redshifted out of resonance with the stationary atoms, hence providing an 
escape route for the photons. The ratio of broad to narrow Ha (and Lya) emission is a function of the shock velocity (HM07; H07); 
it also depends on factors like the pre-shock neutral density and the degree to which the temperatures of the electrons and ions are 
equilibrated. The contribution from the broad Ha line dominates when the shock velocity is < 3000 km s _I and when the narrow 
Ha line assumes Case A conditions (HM07). Existing observations of Ha and Ly/3 emission from 8 Balmer-dominated SNRs are 
catalogued in TableQ] At least 5 of these SNRs are believed to have resulted from Type la explosions. Only SNR 1987A has a clear 
core collapse origin; it is also the youngest SNR in the sample. 

To convert Ha line fluxes to broad Lya luminosities, we use 

^Lya = 47Trf 2 f H « " XT — F Lyo ,/ H ff, (3) 

1 + Aft,, 

where d is the distance to the SNR and %/,„ — 1 is the observed ratio of broad to narrow Ha emission. The quantity r Lya / Ha is the 
ratio of Lya to Ha luminosities (see $3). For SNRs in the LMC, we adopt d — 50 kpc. In the case of the LMC remnant 0509 — 67.5, 
%bn 1S unavailable, so we quote an upper limit. If several values for the Ha flux are given, we simply choose the brightest one (e.g., 
different emission knots of Kepler's SNR). For SNR 1987A, we take the observed value of %tm ~ 1 (Heng et al. 2006), as opposed 
to the theoretically calculated one (~ 0.1; HM07); we use the measured Ha flux to obtain an estimate for the broad Lya luminosity, 
as the measured Lya flux is subjected to resonant scattering. The Cygnus Loop is excluded from our sample due to its low shock 
velocity of ~ 250 km s -1 . In using equation (f3]), we note that the measured Ha and Ly/3 fluxes are mostly from limb-brightened 
portions of the SNRs. Assuming spherical remnants, the intensities from these parts are brightened by factors ~ {R/l a ) 1 ^ 2 (Chevalier 
& Raymond 1978; H07), where R ~ 1 pc is the typical radius of the SNR and l a ~ 10 15 cm is the length scale for atomic interactions 
(assuming density ~ 1 cirT 3 and velocity > 1000 km s _1 ). Hence, the luminosity inferred might be over-estimated by a factor ~ 50. 



1 G07 recorded FUSE spectra only in the 905 — 1 100 and 987 — 1 180 Abands. 
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Fig. 1. Luminosity ratios of Lya to Ha, Lya to LyB, and LyB to Ha, denoted by FLya/Ha, ^Lya/Ly/3 and FLys/Ho-, respectively, as a 
function of the shock velocity, v s . 



3. THE Lya/Ly0 & Lya/Ha Ratios 

The ratio of Lya to Ha luminosities (Fig.[TJ is computed using the methods developed by HM07: 

t{RE,nl + Rt'jiI) + Rt'hI 

T L y a /lia {tll,nl)= — — — , (4) 

e\K-E,n'l' + K T",n'l') + K T' ,n'l> 

where e = Pj /Pj. The quantity Pj is the probability for pre-shock atoms (found in a beam, i.e., at one velocity) to engage in 
charge transfer reactions with ions (thereby creating broad neutrals), while Pi is the probability for the broad neutrals to be ionized 
by both electrons and ions. Physically, broad Ha emission is produced in two ways: charge transfer of the pre-shock atoms to excited 
states of broad neutrals (with a rate coefficient, in cm 3 s , of Rr*,nl)'< creation of broad neutrals in the ground state, followed by 
excitation {Rem) and/or charge transfers between them and ions to excited states (Rt-m)- Hence, e is a measure of how efficient the 
first contribution is relative to the second one. At low shock velocities (v s < 1000 km s _1 ), e > 3 — charge transfer to the ground 
state is the dominant process, and it is efficient to create broad neutrals that subsequently get excited. We emphasize that equation 
is only valid in the case of optically-thin plasmas. 

For Lya, we consider charge transfers (with protons) and excitations (by electrons and protons) to the sub-levels 2p, 3 s and 
3d. For Ha, we consider the same processes, but for the sub-levels 3s, 3p and 3d. Hence, we compute Fi,y a /na(nl,n'l') for nl - 
2p + 3s + 3d and n'l' = 3s + 3d + Bj, lh 2s3p, where the factor of B^ p ,2s = 0.1183 is the fraction of radiative decays from 3p that 
result in Ha, with the remainder going to LyB. For rL ya /iyp(nl, n'l'), we consider instead n'l' = (1 - B-i p ^s)3p. Cascade contributions 
from higher levels are < 5% effects. For example, contributions to Ha from n — 4 are at most ~ (3/4) 3 Z?4 !j 3 p Z?3 / , i 2.s ~ 2%; other 
contributions from 4p, Ad and 4/ are at the < 1% level. 

One can calculate the luminosity ratios for Lyy/Ha, HBfHa and Pa/Ha as well. However, the cross sections for impact excitation 
of hydrogen atoms by protons to the sub-levels 4s, 4p, 4d and 4/ are unavailable at the time of writing^. The cross sections for 



2 We note that Martin (1999) has computed the cross sections for impact excitation of hydrogen atoms by protons to the sub-levels 4s, 4p, 4d 
and 4f, but only for energies of 30 to 200 keV. 
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Fig. 2. Luminosity ratios of Lyy to Ha, H/3 to Ha, and Pa to Ha, denoted by rL yr /Ha, ^Hj3/Ha and rp a /n a , respectively, as a function 
of the shock velocity, v s . Only charge transfers to excited states are considered in these ratios, so they should be used with caution; 
we use only the luminosity ratios for v, > 5000 km s _1 . 



charge transfers to these excited states, however, are available. At v, > 5000 km s _1 , e < 0.5, and we may obtain luminosity ratios 
for Lyy/Ha, H/J/Ha and Pa/Ha to within a factor of 2 (Fig. [2]). A list of the relevant radiative decay fractions, is given in 

Table[3](see Appendix [A] for details). In principle, if the charge transfer and excitation cross sections are known to higher levels in 
the relevant velocity range, one can calculate the luminosity ratios for other lines in the Balmer, Lyman, Paschen and other series of 
hydrogen. 

We use the atomic cross sections of Balanga, Lin & Feautrier (1998), Barnett et al. (1990), Belkfc, Gayet & Salin (1992), Harel, 
Jouin & Pons (1998) and Janev & Smith (1993), as well as those found in the NIST Electron-Impact Cross Section Database. Details 
concerning the cross sections are given in Appendix iBl where we provide fitting functions to them. We consider a pure hydrogen 
gas and include charge transfer, excitation and ionization events between hydrogen atoms, electrons and protons. We employ the 
thin shock approximation, such that the relative velocity between atoms and ions is 3v s /4; this has been shown by H07 to be an 
excellent approximation. At the shock velocities considered, 500 < v s < 10,000 km s , the significance of impact excitation by 
electrons is comparable to that by protons and cannot be neglected. We do not consider broad emission from within the shock front 
(see Appendix ICt . 



4. Results 

The luminosities ratios T^ya/ua, TLycr/Ly/j and rLy/?/Hff are shown in Fig.Q] In the shock velocity range 1000 < v s < 4000 km s~', the 
differences in rL ya /n a and Tiy^/n a between the /3 = 0.1 and 1 cases are factors ~ 2, and are due to the sensitivity to temperature of 
impact excitation and ionization of hydrogen atoms by electrons. This may present a direct and unique opportunity to measure j3. 
We emphasize that our calculations are only valid for the broad lines; the narrow lines have optical depths < r < 1 and Lyman line 
trapping is a non-negligible effect (Ghavamian et al. 2001, 2002). For example, narrow Ly/3 photons may be converted into narrow 
Ha photons and two-photon (2y) continuum. In addition, narrow Lya cannot propagate easily through the interstellar gas. 

We use the data in Table Q] to compute the expected luminosity of Lya, L^ ya (Fig. 0). In estimating a range for LLy«, we only 
consider the observational error bars in Fu a (if available) and allow for a generous range in temperature equilibration between 
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Fig. 3. Left: Expected Lya luminosities, LL ya , from the SNR sample as a function of the shock velocity, v s . Right: Liy a plotted 
versus the age of the SNR. 



electrons and protons, 0.1 < /3 < 1, where /3 = T e /T p . Hence, the displayed error bars for L^yo are not formal ones. We are aware 
of the recent work by Ghavamian, Laming & Rakowski (2007), who showed that there is an empirical correlation between /3 and 
v s — namely, f3 = 1 for v s < 400 km s _1 and f3 oc y~ 2 for v s > 400 km s _1 . For the LMC remnants detected in Ly/3 by G07, we 
compute the range in L\^ a by considering both the Ha and Ly/3 fluxes. We note that the computed (2.96 + 0.05) x 10 35 erg s _1 value 
for broad Lya in SNR 1987A is comparable to the ~ 10 36 erg s figure predicted by Michael et al. (2003). Note that the condition 
Ll/3/Lho > An a / A^p ~ 6.4 is not true in general. This is because the cross section for charge transfers to the level 3p falls below that 
to 3s at a relative velocity ~ 2000 km s _I (Fig. |B.2| i. 

In Table |2j we make some predictions for the Ly/3, Lyy, H/3 and Pa luminosities. It is puzzling that the theoretically expected 
Ly/3 luminosities are about 10 to 20 times higher than those inferred from the observations of G07. In other words, the observed 
Ha fluxes in the LMC SNRs are comparable to the observed Ly/3 ones. We are not certain why this is the case, but we note that 
Ly/3 is more susceptible to absorption by interstellar dust than Ha, and we suspect this effect to play at least some part in the 
discrepancy. Moreover, the Ha and Ly/3 observations were taken at different epochs (Tuohy et al. 1982 versus Ghavamian et al. 
2007). As described in $3j we are only able to provide rough predictions for Lyy, H/3 and Pa, and only in the cases of 0509 — 67.5 
and SNR 1987 A, as these SNRs have shock velocities > 5000 km s~' . 

We can make some estimates for the expected 2y continuum as well, which is produced in the 2s — > Is transition. In the case 
of an optically-thin plasma, the 2s — > 2p transition is negligible as collisions are unimportant. In Table [2] we make conservative 
predictions for the 2y continuum luminosity from both broad and narrow atoms, but consider only charge transfers and excitations 
to the 2s level. Additional contributions from n — 3 range from ~ (2/3) 3 B3 Pi 2.s ~ 4% (Case A) to ~ (2/3) 3 ~ 30% (Case B); those 
from n = 4 are «: 1%. We only wish to make the point that L-i y is comparable to Li^ a and Lu a , and thus the 2y transitions are a 
potentially observable source of continuum. In the case of Galactic and LMC SNRs, the Galaxy Evolution Explorer ( GALEX) is in 
principle able to measure the low-frequency wing of 2y decay, using its 135 — 175 and 175 — 280 nm channels. By comparing Ha 
and 2y emission, it will be possible to directly estimate emission from the SNR shock due to broad Lya (and the contribution of 
narrow Lya that cannot reach us). This is an additional, unique source of information on the detailed physical processes in shocks. 

Several sources of uncertainty can affect the predicted luminosities. These include uncertainties in the age of the SNR, f age , the 
distance to it, d, the measured non-radiative component of the Ha flux, the temperature equilibration between electrons and ions, 
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Fig. 4. Luminosity ratios of the 2y continuum to Ha, r2 7 /n a - Only charge transfers and excitations to the 2s level are considered 
(see text). 



and the atomic cross sections used. Uncertainties in the cross sections are typically ~ 10%. For charge transfer to excited states, the 
uncertainty can be as much as 30% (R.K. Janev 2007, private communication). The predicted luminosities have not been corrected 
for reddening by dust. 



5. Discussion 

SNR 1987A is a unique example of a B aimer-dominated SNR. By virtue of adiabatic expansion cooling, the SN ejecta comprises 
mostly neutral hydrogen; it rushes out at velocities > 12, 000 km s _1 (Michael et al. 2003; Heng et al. 2006). The non-radiative Ha 
and Lya result from the interaction of the ejecta with the reverse shock and not the blast wave (Heng 2007). As SNR 1987A has a 
Type II origin, it is possible to produce Balmer and Lyman lines via this mechanism; this is obviously not possible with Type la's. 
Smith et al. (2005) have predicted that the Ha and Lya emission from the reverse shock of SNR 1987A is shortlived (~ 2012 to 
2014) and will be extinguished by the increasing flux of extreme ultraviolet (EUV) and X-ray photons traveling into the pre-shock 
region and ionizing the atoms — pre-ionization. This is marginal evidence that broad Lya from SNRs of a core collapse origin will 
be short-lived, i.e., < 100 years. In general, for this scenario to work, some interaction of the blast wave with the ambient material 
is needed, but if it is too strong the pre-shock gas becomes ionized (R. Chevalier 2007, private communication). 

To further investigate the viability of the short-lived, non-radiative Lya hypothesis, we examine the sample of optically identified 
SNRs by Matonick & Fesen (1997), who studied an ensemble of 12 SNR samples from different galaxies, including the Small 
Magellanic Cloud (SMC), LMC, M31 and M33, with distances up to 7 Mpc. In galaxies like NGC 2403, M81 and M101, the SNRs 
are associated with star-forming regions and most of them probably have a Type Ib/c origin. In most cases, the measured Ha flux is 
~ 10~ 15 erg cirT 2 s _1 and the inferred luminosity is ~ 10 36 erg s _1 . Since Matonick & Fesen (1997) did not provide Ha line profiles, 
it is impossible to estimate the proportion of the Ha emission that is non-radiative. Furthermore, their selection criterion is based on 
picking out objects with [S n]/Ha > 0.45, which will not detect SNRs with predominantly non-radiative Ha emission. 

Shull & Silk (1979) computed the temporally-averaged Lya luminosity from radiative shocks of a population of Type II SNRs, 
assuming low metallicities, to be 

L SS79 = 3 x 10 43 erg s-^V^sn, (5) 
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where TVsn is the number of supernovae (SNe) a year. They considered SNRs in both the ST and the PDS stages, and v s = 20 to 120 
km s _1 . Chariot & Fall (1993) remark that the numerical coefficient in the preceding equation is about 40% lower if one assumes 
solar metallicity. 

A very conservative upper limit on the broad Lya from the Matonick & Fesen (1997) samples can be obtained if one generously 
allows for all of the Ha to be broad, for the shock velocities to be low (~ 500 km s _1 ) such that F LyQ ,/ H a ~ 100, and for the 
non-radiative emission to last ~ 10 4 years. Even in this very unlikely scenario, L Lya ~ 10 42 erg s _1 is only about 0.1Lss79- Hence, 
our charge transfer mechanism is not energetically competitive. There is the possibility a SNR can produce both radiative and 
non-radiative components of Ha. Well-known examples are Kepler (Fesen et al. 1989; Blair, Long & Vancura 1991) and RCW 
86 (Long & Blair 1990; Smith 1997). There is also the possibility that the non-radiative emission from the SNR is inhibited. For 
example, Foster (2005) observed and studied the Galactic SNR 3C 434.1 (f age « 25, 000 yr; d = 4.5 + 0.9 kpc; possible Type Ib/c), 
which formed inside the eastern portion of a pre-existing stellar-wind bubble of interior density ~ 0.1 cm 3 . Strong Ha emission 
(6.1 + 0.4 x 10 36 erg s _1 ) is measured from the eastern side; it is believed to be from a radiative shock. Being farther away from 
the western wall of the bubble, the shock on the western side is essentially still in free expansion and produces no measurable, 
non-radiative Ha. 

Our SNR sample and the considerations of SNR 1987 A lead us to believe that if the short-lived emission contribution from Type 
Ib/c and Type II SNRs in young galaxies exists, it has a luminosity of 

Lhya,CC ~ 10 38 erg S _1 f e mit,2#SN, (6) 

where f em i t = f e mit,2 1 00 years is the length of time we expect core collapse SNRs to produce shock-induced Lya emission. On the 
other hand, thermonuclear SNRs are expected to have f em i t = / e mit,4l0 4 years ~ f PD s. However, they are also believed to be much 
scarcer at high redshifts. For example, Dahlen et al. (2004) estimate that only 5% to 7% of available progenitors explode as Type la 
SNRs. Therefore, the expected luminosity is 

LhyaM ~ 10 38 erg S _1 W,4#SN,-2, (7) 

where Nsn-2 is the number of SN per year in units of 0.01. We conclude that for both core collapse and thermonuclear SNRs, 
the expected luminosity from broad Lya is only a ~ 0.001% effect, compared to the mechanism of Shull & Silk (1979). Lya line 
luminosities from z ~ 3 to 5 galaxies have been observationally determined to be ~ 10 42 to 10 43 erg s _1 (e.g., Saito et al. 2007), in 
general agreement with theoretical expectations. In addition, the lifetime of an emitting atom is approximately the length of time 
corresponding to one atomic length scale, and is only f m f p ~ l a /v s ~ 10 7 n 1 v^g s, where v. Sj 8 = v s /1000 km s _1 (H07). 

We have restricted our analysis to homogeneous and static media. Though broad, non-thermal Lya emission has never been 
observed, these photons are produced in SNRs and hence the non-radiative Lya luminosity is a part of the intrinsic Lya spectrum 
of young galaxies. The optical depth for a broad photon in the line wings is (Verhamme, Schaerer & Maselli 2006) 

t ~ 0.26r 4 1/2 ^V H; 20V W >12.85, (8) 

where T = 10 4 T 4 K and N H = N H ,2ol0 2() cm 2 are the temperature and obscuring hydrogen column density of the medium, 
respectively. The turbulent velocity in the interstellar medium is b = 12.85^12.85 km s _1 (Verhamme, Schaerer & Maselli 2006), 
while v w = lOOOv^g km s _1 is the velocity of the emitting atom in the line wings. Multiple scattering is important for t > 0.3 
(Chevalier 1986) and any realistic treatment of non-thermal Lya lines in a young galaxy has to include radiative transfer effects, 
which we have neglected in our analysis. 
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Table 1. Ha & Ly/3 Observations of SNRs 



Object 




d 




2vi>n 


Ho- 1 


LyyST 


SN Type 


References 




(yr) 


(kpc) 


(km s- 1 ) 




(erg s" 1 ) 


(ergs ') 






0505—67.9 


~ 10,000 


50 


464 — 744 


> 0.7 


7.15 34 


(4.58 ± 0.23)34 


la 


5,7,10,12 


0509—67.5 


> 250 


50 


5200—6300 




< 3.29 34 


< (2.33 ±0.18)34 


la 


5,7,9,10,12 


0519—69.0 


>400 


50 


1032—1809 


0.8 ±0.2 


3.19 34 


(2.93 ± 0.07)34 


la 


5,7,9,10,12 


SN 1006 


992 


2 +0 - 4 


2290—3111 


0.84^ 


(4.20 ± 0.92)3i 




la 


4,7 


Kepler 


384 


2.9 ±0.4 


1518—2446 


0.72 ± 0.37 


6.74 31 






1,7 


RCW 86 


1802* 


2.5 


496—662 


1.18 ±0.03 


(2.09 ± 0.25)30 






3,7,9,13 


SNR 1987A 


18 


50 


7840—9200 


1? 


(2.98 ± 0.33)34 




11 


6,7,11 


Tycho 


406 


1.5—3.1 


1631—2344 


0.67 ±0.1 


(4.15 ±2.31)30 




la 


2,3,7,10 



f Inferred broad line luminosities from published line fluxes. Note that At is shorthand notation for Ax 10 fc . 
t Assuming that RCW 86 is the remnant of SN 185. 

* Quoted for 0.1 < P < 1, where = T e /T„ is the ratio of electron to proton temperature. 

1: Blair, Long & Vancura (1991); 2: Chevalier, Kirshner & Raymond (1980); 3: Ghavamian et al. (2001); 4: Ghavamian et al. (2002); 5: Ghavamian et al. (2007); 6: Heng et al. (2006); 7: 
Heng & McCray (2007); 8: Long & Blair (1990); 9: Rest et al. (2005); 10: Smith et al. (1991); 1 1 : Smith et al. (2005); 12: Tuohy et al. (1982); 13: Zombeck (1982). 



Table 2. Predicted broad Lye, Ly/3, Lyy, H/3, Pa and broad/narrow two-photon luminosities (erg s ) 



Object 


Lya 


Ly/3 


Lyr 






Pa 


Broad 2y 


Narrow 2y 


0505—67.9 


(4.90 


±4.61)36 


(1.09 


± 0.34)36 












(6.14 


± 1-38)35 


(2.51 


±0.12)35 


0509—67.5 


(2.68 


± 2.05)35 


(1.29 


±0.14) 3 5 


(2.04 ± 0.08)34 


(5.97 


± 0.36)33 


(1.07: 


t 0.06)33 


(6.16 


± 0.02)35 






0519—69.0 


(1.46 


± 1-29)36 


(2.93 


± 1.45)35 












(5.00 


± 1-37)35 


(2.86 


± 1-83)35 


SN 1006 


(8.42 


± 2.57)32 


(2.54 


±0.15)3? 












(5.91 


± 2.04)32 


(2.97 


± 2.03)32 


Kepler 


(2.95 


± 1-24)33 


(5.30 


±2.18)32 












(1.01 


±0.15)33 


(5.82 


± 3.70) 3 o 


RCW 86 


(2.30 


± 0.94) 12 


(3.48 


± 1.47)3i 












(1.76 


±0.48)3i 


(4.11 


± 0.06)30 


SNR 1987A 


(2.96 


± 0.05)35 


(7.64 


± 0.59)34 


(1.99 ±0.44)34 


(5.96 


± 1-40)33 


(1.06: 


t 0.25)33 


(6.43 


± 0.52)35 


(5.75 


± 0.70)35 


Tycho 


(1.03 


± 0.04) 3 , 


(2.32 


±0.44)n 












(6.83 


±4.52)3, 


(5.38 


±4.78)3, 



Note: Af, is shorthand notation for A X 10 /j . 
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Table 3. Radiative decay fractions 



Quantity 


Value 


Relevance 


-83/7,2 j 


0.1183 


Ha, Ly/3 




0.5841 


H/3, Pa 


-64.1,3/) 


0.4159 


H/3 Pa 


Ba p ,\ s 


0.8402 


H/3, Lyy, Pa 




0.1191 


H/3, Lyy, Pa 


#4p,3s 


3.643_ 2 


H/3, Lyy, Pa 


Bnp^d 


4.282^3 


H/3, Lyy, Pa 


B\d,2p 


0.7456 


H/3, Pa 


Bn,d,ip 


0.2544 


H/3, Pa 



Note that A b is shorthand notation for A X l() b . 



Appendix A: Ratio of Einstein A-coefficients 



To compute the rate coefficients for Lyes', Ly/3, Lyy, Ha, HB and Pa, one needs to calculate the ratio of Einstein A-coefficients. The 
Einstein A-coefficient for hydrogen, A^n'i't is the radiative decay rate (s -1 ) from the levels nl to n'V . The radiative decay fraction is 



B„. 



Ln'l' 



^ A„l,n"l" 



\n"l" 



(A.l) 



where the sum is over all transitions nl — > n"l" permitted by the electric dipole selection rule, /" = l±\. A list of the relevant radiative 
decay fractions is listed in Table|3] For example, to compute Tiyj/Hainl, n'V), we need to consider nl - B^^p^s + B^p^Ap + B^xAd- 
Our computed value of Bi p 2 S = 0.1 183 is in close agreement with the 0. 1 184 value quoted by Martin (1999). 

The Einstein A-coefficients are proportional to the square of the magnitude of the radial integrals, \R"J \ 2 . (See Appendix A2 of 
HM07 for details on how to calculate them analytically.) As a check, we have compared our computed values of \R n J | 2 to the ones 
tabulated by Green, Rush and Chandler (1957), and find them to be in agreement. 



Appendix B: Atomic Cross Sections 

Cross sections for interactions between protons and hydrogen atoms (charge transfer and excitation) to the sub-levels 3 s, 3 p and 
3d are computed in Balanfa, Lin & Feautrier (1998) and kindly provided to us by C. Balanga (2007, private communication). We 
note that these calculations utilize a two-center atomic-orbital (TCAO) close-coupling method with an asymmetric (TCAO-A) basis 
set of 26 states, so as to avoid spurious numerical oscillations caused by using a traditional, symmetric set (TCAO-S). We fit these 
cross sections, as well as those from BeMc, Gayet & Salin (1992), Harel, Jouin & Pons (1998) and the NIST Electron-Impact Cross 
Section Database using the function: 



T (x; A) = exp 



8 \ 



= 1 



(B.l) 



where the coefficients A = A,- for < i < 8 are the fitting parameters. The quantities C, are the Chebyshev orthogonal polynomials: 



C\ (x) = x, 


(B.2) 


C 2 (x) = 2x 2 -1, 


(B.3) 


C 3 (x) = 4x 3 - 3x, 


(B.4) 


C 4 (x) = 8(x 4 -x 2 )+l, 


(B.5) 


C 5 (x) = 16x 5 - 20x 3 + 5x, 


(B.6) 


C b (x) = 32x b - 48x 4 + 18jc 2 - 1, 


(B.7) 


d (x) = 64x 7 - 1 12x 5 + 56x 3 - Ix, 


(B.8) 


Cg (x) = 128x 8 - 256x 6 + 160x 4 - 32x 2 + 1. 


(B.9) 




(B.10) 



The fitting variable x is defined as 
In (E/E min ) - In (E m:ix /E) 



(B.ll) 



In (Emax/ E m [n) 

where E is the relative energy between the collidants; E m i n and E m . dx are the respective minimum and maximum energies to which the 
data are available. We use the Levenberg-Marquardt algorithm, which combines the steepest descent and inverse-Hessian function 
fitting methods, as implemented in IDL. The fits are sensitive to the initial values of the parameters fed to the algorithm; we use the 
values of the fit parameters for cr T p 2 S in Barnett et al. (1990) as a guide. In providing "measurement errors" to our fitting algorithm, 
we assume a fiducial error of 10%. Selected cross sections are shown in Figs. lB. lllB.2IIB.3llB.4l and lB.5l while the fitting coefficients 
are presented in Table IB. ll 
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Fig. B.l. Cross sections for charge transfers between hydrogen atoms and protons, to the sub-levels 2s and 2p, taken from Belkfc, 
Gayet & Salin (1992) and Harel, Jouin & Pons (1998). Also shown are the fitting functions of Barnett et al. (1990), denoted by 
"B90". 



Appendix C: Broad Emission from Within the Shock Front? 

We have considered only the case of hydrogen atoms crossing the shock front (in the frame of the front) and interacting with 
protons. In this sub-section, we examine the possibility of Lya being created within the shock front. The width of the shock front 
in collisional shocks is on the order of an atomic mean free path, l m f p , assuming a pure hydrogen gas. Zel'dovich & Raizer (1966) 
have shown that for weak shocks, the collisional shock width is 

<5~/mfb — — — , (C.l) 
pi -po 

where po and p\ are the pre- and post-shock fluid pressures, both in the case of a viscous shock with no heat conduction and for a 
heat-conducting shock with no viscosity. When the change in pressure across the shock front is comparable to the magnitude of the 
pre-shock pressure, 6 ~ l m f p (see also Landau & Lifshitz 1963). Even in the limit of infinite Mach number, Sakurai (1957) finds that 
tf/Znrfj, = 1.42. 

The question is: how robust is the assumption of shocks in B aimer-dominated SNRs being collisionless? This occurs when the 
electron and proton gyroradii — r e and r p , respectively — are much smaller than Z mtp . The typical value of the magnetic fields in 
SNRs is B = fi- 4 l(r 4 G. For example, Strom & Duin (1973) find 3 x 1(T 4 and 5 x 1(T 4 G for Tycho and Cas A, respectively. The 
electron gyroradius is 

r e ~ 10 5 cm v e ^B~_\, (C.2) 

where v e $ = v/1000 km s" 1 is the velocity of the electron. For protons, we have r p ~ 10 8 cm v p $/B-4. It follows that the transition 
from collisionless to collisional shock occurs when the gyroradii ~ / m f p , or when the density of particles is n e > W W v~\cr a \ L 15 #-4 
crrT 3 and n p > 10 7 v p 1 g cr^ 1 15 B_4 crrT 3 , where cr a = cr fl _ 15 1CT 15 crrT 2 is the typical value of the cross section for atomic interactions 
(charge transfer and ionization; see HM07 and H07). These densities are much larger than typical values for the interstellar medium 
(~ 1 cm" 3 ) or even for molecular clouds (~ 100 to 1000 cm" 3 ). 
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Fig. B.2. Cross sections for charge transfers between hydrogen atoms and protons, to the sub-levels 3 s, 3 p and 3d, taken from 
Belkfc, Gayet & Salin (1992) and Harel, Jouin & Pons (1998). Shown for comparison are the calculations of Balanca, Lin & 
Feautrier (1998), denoted by "BLF98". 



About exp (— /j//mfp) of the hydrogen atoms cross the shock front without being ionized, where l s is the width of the shock. In 
collisional shocks, l s ~ 6. In collisionless shocks, l s ~ r e , l s /l m f p <K 1 and virtually all of the atoms pass through. We thus conclude 
that broad Lya is probably not produced in a significant amount within the shock front, consistent with the findings of H07. 
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Fig. B.3. Cross sections for charge transfers between hydrogen atoms and protons, to the sub-levels 4s, 4p, Ad and 4f, taken from 
Belkfc, Gayet & Salin (1992) and Harel, Jouin & Pons (1998). 

Table B.l. Fitting Parameters for Various Cross Sections 



Reaction 


A 


A, 


A 2 


A 3 


A 4 


A 5 


A 6 


A 7 


A 8 


°" T,p,2s 


-9.76244, 


-8.02645 


-9.03528 


-6.82148-, 


1.41138 


-8.43613_ 3 


-1.31162-1 


-3.07135-1 


3.60166-1 


°" T,p,2p 


-9.47856, 


-1.36982, 


-7.48411 


-6.55633_i 


1.08973 


3.45255_i 


-2.02178.! 


-2.99212-1 


2.88704_i 


°"T,p,3s 


-1.03831 2 


-5.12682 


-1.03311, 


-2.36941_i 


7.99974_i 


6.55927_i 


-2.64604-1 


-6.82557_i 


7.65819-1 


°" T,p,3p 


-1.04117 2 


-9.31370 


-1.09011, 


2.72437_i 


8.33272-1 


9.72800-1 


-3.40175-1 


-8.21497.2 


1.06280-1 


f T,p,3d 


-1.06187 2 


-8.82542 


-1.05048, 


9.29637_, 


9.65598_i 


7.99638.3 


1.21539-1 


-6.52001_i 


4.15571-1 


f T,p,4s 


-1.050182 


-5.69159 


-9.63683 


-8.65652_i 


1.06090 


9.33260-1 


-8.64227_i 


-1.30903-1 


4.32513-1 


C T,p,4p 


-1.06245 2 


-7.18909 


-1.07639, 


1.84045-1 


1.10514 


3.64508., 


-2.30906-1 


-4.14293.1 


4.42254_i 


C T,p,4d 


-1.08456 2 


-7.54281 


-1.11995, 


1.15545 


9.55810-1 


9.05249.2 


-1.37453., 


-4.09046.1 


3.40271-1 


f T,p,4f 


-1.10117 2 


-8.87472 


-1.00318, 


3.41036_i 


1.43444 


1.22887-2 


-4.35737-1 


-7.29281.2 


1.87428., 


f E,e,2p 


-7.67030, 


-9.47539., 


-6.80475_, 


2.77805-1 


-1.18599., 


6.65029-2 


-4.38768_2 


2.20102- 2 


-2.24010.2 


C E,e,3p 


-8.02681, 


-1.00579 


-6.24421_, 


2.59918-1 


-1.13262., 


6.48909-2 


-4.34274.2 


2.15313-2 


-2.31916.2 


C E,p,3s 


-8.43473, 


2.69180 


-8.94871_, 


-1.25347-1 


-1.27427-, 


2.16335-, 


7.34001-3 


-2.19545-1 


-1.76763., 


C E,p,3p 


-8.10541, 


2.27478 


-4.95559_, 


1.56523_ 2 


2.45763.2 


-1.29970-1 


-2.66319-2 


-6.13906-2 


2.46815-, 


C E,p,3d 


-8.22914, 


7.43945-, 


-5.04030-, 


5.76290-2 


8.40403-2 


-3.04995-1 


7.62616-2 


2.56426-1 


-3.65589-1 



Note: A b is shorthand for A X 10 6 . The following subscripts are used: "T" (charge transfer), "E" (impact excitation), "p" (proton) and "e" (electron). 
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Fig. B.4. Cross sections for impact excitation of hydrogen atoms by protons, from the two-center atomic-orbital (TCAO) close- 
coupling calculations of Balanca, Lin & Feautrier (1998). Shown are the fits to the TCAO-A calculations. 
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Fig. B.5. Cross sections for impact excitation of hydrogen atoms by electrons, to the sub-levels 2p and 3p, from the NIST Electron- 
Impact Cross Section Database. 



